Fermi-like behavior of weakly vibrated granular matter 
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Abstract 

Vertical movement of zirconia-yttria stabilized 2 mm balls is measured by a laser facility at the surface of a vibrated 3D 
granular matter under gravity. Realizations z(t) are measured from the top of the container by tuning the fluidized gap with 
a ID measurement window in the direction of the gravity. The statistics obeys a Fermi-like configurational approach which 
is tested by the relation between the dispersions in amplitude and velocity. We introduce a generalized equipartition law to 
characterize the ensemble of particles which cannot be described in terms of a Brownian motion. The relation between global 
granular temperature and the external excitation frequency is established. 



A vibrated granular medium (GM) exhibits a wealth of 
intriguing physical properties 0, 0, EJ 0, @ • Since energy 
is constantly being added to the system a nonequilibrium 
steady state (s.s.) is reached In we have studied 
the spectrum properties of vibrated GM under gravity, 
and shown that in the weakly excited regime the dynam- 
ics of the fluidized particles cannot be described as simple 
Brownian particles, this fact leads us to the conclusion 
that in order to describe the cooperative dissipative dy- 
namics of the GM particles, it must be done in terms of 
generalized Langevin particles 0, . 

Recently Hayakawa and Hong [Jjj introduced the ap- 
proach of thermodynamics of a weakly excited granu- 
lar matter, in particular vibrated GM by mapping the 
nonequilibrium system with a "Fermion like" theory. Our 
experimental conditions allow us to consider a N parti- 
cles system of n-rows in a cylindrical container as a ID 
degenerate Fermi system. Two experiments based on a 
laser were set up to investigate the occupation dynam- 
ics at the fluidized gap of the n-rows GM. The first one 
considers a realization z(t) of one particle from the top 
of the container by tuning the fluidized gap with a ID 
window in the gravity direction, see Fig. 1(a). It is clear 
that these realizations mainly correspond to macroscopic 
Fermi-like particles (mFp) from near the "Fermi level". 
The second one concerns the measurement after a long 
integration time of the mass profile which corresponds to 
the Fermi-like profile, see Fig. 1(b). By focusing on the 
configurational properties of an excluded volume theory, 
the s.s. mass profile can be understood in terms of a con- 
figurational maximum principle assumption. Excluded 
volume interactions of the GM do not allow two grains 
to occupy the same state (gravitational energy) , thus the 
number of configurations is W = Y[i [f2!/ATj!(f2 — Ni)\]. 
Following Landau, to study a non-equilibrium system, 
the maximization of S = InW yields that the profile 
is </>(e) = [1 + Q exp(/3e)] 1 £1/ mgD where m, D are the 
mass and diameter of our balls; is a Lagrange multiplier 
parameter, Q^ 1 = exp(/3/x) the fugacity, and e = mgDs 
with s = 0, 1, 2, 3 • • • . Introducing the normalization: J 
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FIG. 1: Distance measurement (amplitude realizations z(t)) 
corresponding to 13 layer GM (h — 21mm) with 1.99mm 
Z r O% — Y2O3 balls in a glass container of 30mm diameter. 
Set up for (a) the single particle measurement and (b) the 
laser light barrier experiment. 



(j)(e)de = N, we get exp(N(3mgD/n) = 1 + exp(/?/x). 
Therefore the zero-point "chemical potential" is fio — 
mgDN/Cl, where N/£l is the number of balls in an el- 
ementary column of diameter D. From these considera- 
tions it is trivial to see that without vibration the centre 
of mass (cm.) is characterized by z c m = fi /2mg = h/2. 
For vibrated GM the Lagrange parameter (3 is a non- 
trivial function of the velocity fluctuations. Before going 
ahead, let us denote <f>(e)/ N as the cumulative probabil- 
ity that the energy e will be occupied in an ideal GM layer 
at the nonequilibrium s.s. characterized by the global 
temperature (3~ l . 

From cf)(e) it is possible to calculate the cm. expan- 
sion as a function of f3, the mean energy per particle, its 
square dispersion of, etc. Many questions concerning the 
GM layer system are still open, in fact the stochastic mo- 
tion of the fluidized particles is not entirely known |ll| . 
For example it is important to test that the spectrum of 
the realization S z (f) does not behave as Brownian par- 
ticles (l// 2 ), but it has a more complex behavior, 1//", 
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related to a cooperative dynamics 0,0 ■ Thus an exhaus- 
tive analysis of the realizations z(t) of these macroscopic 
Fermi-like particles should be made. We have measured 
z(t) and we calculate the Lagrange parameter [3 to show 
its complex relation to the kinetic energy. 

Amplitude dispersion vs. the velocity dispersion 
A sinusoidal vibration is driven by a vibration plate on 
the GM bed (with intensity T — Auj 2 /g, where A is the 
amplitude, g is the acceleration of gravity and lo = 2?:f e 
the frequency of the plate). The vibration apparatus 
is set up by an electromagnetic shaker (TIRAVIB 5212) 
which allows 0] for feedback through a piezoelectric 
accelerometer the control of / e and T in the range of 
10- 7000Hz, and 2 — AOg respectively. The control loop 
is completed by an Oscillator Lab-works SC121 and a 
TIRA19/z amplifier of lkw. The n-rows GM bed set 
up were Z r Oi — Y2O3 balls with D = 1.99mm and 
m = 26.8 ± O.lmgr, into a glass container of 30 mm of 
diameter with steel bottom, see Fig. 1(a). The experi- 
ments were carried out in a chamber at latm of air with 
5.8 ± 0.2gr/m 3 of water vapor. The absolute humidity 
was controlled by using a peltier condenser and a control 
loop through a thermo-hygrometer. The humidity is of 
major relevance in order to control the particle-particle 
and particle- wall contact forces 0, El ■ Under such 
humidity controlled conditions, no surface convection or 
convection rolls were observed in the GM, nor rotational 
movement of the bed with respect to the container, which 
is typical for a content of water vapor > 10gr/m 3 . 

The z(t) of one particle was followed in a window 
of 12mm with a laser device by using a triangulation 

method, see Fig. 1(a). A laser emitter with a spot of 
70/mi and a linear image sensor (CCD-like array) enables 
a high speed measurement with 100/isec sampling. The 
linear image sensing method measures the peak position 
values for the light spots and suppresses the perturbation 
of secondary peaks, which makes possible a resolution of 
1/im. The shaker and the laser displacement sensor were 
placed on vibration-isolated tables to isolate them from 
the external vibrations, and the displacement sensor from 
the experiment vibrations. The z(t) is a measure of the 
variations of the distance (difference) between the par- 
ticle and the sensor around the surface of the GM bed 
(fluidized gap). The measured z(t) without excitation 
reveals a white noise < 10/im. Then our set-up effective 
resolution is no higher than 10/im. We have shown0 
that depending on the external excitation the z(t) can 
show from quasi non-erratic parabolas, for the movement 
under gravity, to realizations of larger rugosity. 

The registers of z{t) were taken with a 9354 C Le 
Croy Oscilloscope of 500MHz. The velocity V(t) = 
dzj dt of the mFp was calculated numerically for At = 
100/isec from z(t) registers. The dispersions o~ z = 

^J(z(t) 2 ) - (z(t)) 2 and a\ = (V(t) 2 ) - (V(t)) 2 were ob- 
tained from a window of 2 second for each pair of registers 



{z(t), V(t)}. In Fig. 2(a) we report a z against o v for fixed 
T = 10, 20 and several f e from 60-180Hz for GM beds of 
h = 21mm and 12mm. 

For weakly excited GM the displacement of the flu- 
idized particles, in the gap, can be studied from the pro- 
file <^(e). In fact a nonequilibrium s.s. density P(e = 
mgz), characterizing the motion of the fluidized mFp, is 
sustained by the input of energy from the plate colliding 
periodically with the GM bed; i.e., a current of particles 
near /io — which is proportional to a gradient of </>(e) — will 
be balanced by the random input of mass coming from 
the periodic movement of the plate. It is clear that P(e) 
will be a narrow density around /io, so we characterize 
the movement of the mFp at the Fermi-like sea by 



#(e) , 
P(e) oc — , where e 



mgz = mgDs. (1) 



At the nonequilibrium s.s. the dispersion a z can be calcu- 
lated from P(e = mgz), but a rather simple and analyti- 
cal expression for a characteristic length scale z* can be 
obtained by solving e* from the following consideration 



qP(fio) = P(fio + e*); 0<q< 1, 



(2) 



where q is a "cumulant" parameter. If P(e) were Gaus- 
sian the value q = 1/ yfe would give the exact dispersion 
e* = +a e . We have tested that our conclusions are not 
changed for values q ~ \j\fe.. Using <p{e) in the expres- 
sion for P(e) we get for the characteristic scale e* 
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(2,4-1) 
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2A^q 



(3) 



where A = e^ Mo , then by putting q ~ 1/ i/e in © it gives 
the amplitude dispersion z* = e*/mg as a function of 
P. Now the task is to determine (3 as a function of the 
kinetic energy of the mFp in the fluidized gap. 

If collisions were elastic, in a ID ideal gas the equipar- 
tition theorem says that total kinetic energy per particle 
is related to the Lagrange parameter by may /2 = /3 _1 /2. 
Our conjecture for a weakly excited GM is to generalize 
the equipartition law to 
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(4) 



where m* — 8 m accounts for inelastic factors, and 
AN/ N is a relative factor that counts the thermody- 
namically "active" mFp in the fluidized gap. In fact, a 
granular gas since its non-Gaussian velocity distribution 
reveals an inelastic gas heated in a non-uniform way, with 
the expected high energy tail e - constant v3/2 . 
The factor AN/ N can be calculated from 



AN = 



= 4— 



(e) de = N ( In [2e p ^ - l] - 1 



(5) 
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So the implicit equation to solve (3 is 



K^" ,[2 ' 



-1-1 



= SlTMJy. 



(0) 



Note that in the high temperature limit [3[io <C 1 and 
for the elastic case (5=1 we recover the equipartition 
theorem. This situation is just what we have found ex- 
perimentally for one steel ball in a narrow glass cylinder 
pj. In that experiment, when we compared the relation 
a z vs. ay, we reduced the dissipation during the vi- 
bration and assured that there is no rotation of the ball 
during its movement z(t), then from energetic considera- 
tions: mga z = ^ma Vl predicting a line with slope l/2g. 
The opposite situation is in the limit [3^iq 3> 1, in this 
case we arrive at the Low Temperature (LT) scaling. 



(3- 



a v 
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Due to the fact that dissipation and degrees of freedom 
are functions of the external parameters, we expect that 
the analysis of the complex behavior of vibrated GM will 
be enlightened from the study of a z — a z (oy) ■ Thus 
an important point would be to test experimentally our 
theoretical predictions. Noting that z* = a z = e*/mg it 
is simple to see that at LT gives 



(3-2 



ln(l + VT— D-hiy/q 



mga z 



then using (JJJ) we arrive to 



(ln(l + v ^^)-ln^) 1 



/ 4hS 
~g~W2 



(8) 



(9) 



Thus we got an explicit LT formula a z — a z (cry) as 
a function of the dissipative parameter S, which in fact 
is a function of the external parameters T and f e . In 
Fig. 2 (a) we report the measurement of a z against a v 
for two experimental studies of a GM bed with h = 21 
and h — 12 mm at V = 10 and r = 20 respectively, for 
several f e from 60-180Hz. In that figure we also show 
the fit with our theoretical prediction (resolved per least 
squares) showing a very good agreement for 8 ~ 0.0009 
and S ~ 0.0064 (T = 10 and T = 20 respectively). In 
Fig. 2(b) we show the corresponding (3 against er z , where 
the {(3} data set were calculated from the {cry} experi- 
mental data set for the two experimental studies, using 
(Q). By considering the mass of the ZrO^ — YiO-s ball and 
q — we represent in Fig. 2(b) the log-log plot of the 

equation (jSJ , showing an excellent agreement between the 
experimental data and our theory. The two experimen- 
tal data sets are on the same curve due to the fact that 
we use for the two experiments the same Z r Oi — YiOj, 
balls. Note that if all stochastic realizations could be 
understood in terms of a Brownian oscillating movement 
around /zq, the P(z) would correspond to exp (—z 2 C(3), 
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FIG. 2: (a) a z vs. a v , and the corresponding fittings by using 
e* vs. Oy solved from @, or from ||HJ in the LT approxima- 
tion. The excitation frequency was from 60 — 180Hz. (b) 
Lagrange parameter f3 vs. a z . 



with C a constant and /3 _1 proportional to the temper- 
ature. Then we would have obtained a z cx f3~~ 1 / 2 which 
is not the case reported experimentally in Fig. 2(b). Here 
we point out that in order to describe the spectrum of 
the fluidized particles we should use a non-Markovian de- 
scription 0. Unfortunately we still do not have a time- 
dependent statistical description for the mFp. Global 
temperature against the external excitation. Eq. (JTJ is 
the LT approximation of our generalized equipartition 
theorem for a GM experiment out of equilibrium. Now 
we would like to find a relation between the Lagrange 
parameter [3 and the external parameters characterizing 
the input of energy. The maximum kinetic energy, per 
particle, transferred by the oscillating plate must be pro- 
portional to the effective mass m* and the dimcnsionless 
velocity Auj / \fgd, on the other hand the maximum po- 
tential energy related to the fluidized gap is proportional 
to the variation of the cm. at the global temperature 
Thus we get the relation 



Alo 



mgAz c 



(10) 

iga J no 
Where Az c . m . = U/N - /i /2, with U = J Q °° e<j){e) de. At 
LT we get Az c . m . ~ {it/ f3) 2 j 6/io, then 



0- 



mh\jZgSjd 



Au, h = no/rug 



(11) 



Using the relation Auj — Tg / 2irf e we can transform 111|) 
in terms of the variables that we have fixed in our exper- 
iment. From these considerations it is trivial to see that 
for a given intensity T and increasing frequency f e —fOo 
the "temperature"^ 0, (i.e.,cr 2 — * and z c . m , — > h/2). 
In Fig. 3(a) we present the behavior of the global temper- 
ature as a function of the excitation frequency show- 
ing an agreement with our theoretical prediction. In 
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FIG. 3: Comparison of the two experiments for a GM bed of 
h — 21mm under an acceleration F — 10. (a) the Lagrange pa- 
rameter /3(/ e ) by measuring the realizations z(t) from the top 
of the container, (b) profile <f>(z)/ N integrated from the laser 
light barrier during 3h corresponding to f3 = 4.4 ± 0.2/^J -1 
for f e = 100Hz, and (1 = 10.5 ± 0.5^J _1 for f e = 125Hz. 



Fig. 3 (a) a least squares fitting is also shown giving a 
slope of 0.06/4J -1 s, while dHJ gives O.O^J^s. The 
dispersion in Fig. 3(a) is mainly introduced by the nu- 
merical calculation of dz/dt. The nonequilibrium s.s. 
Fermi-like profile. We have remarked that the profile 
<f)(e) j N gives the probability that the energy e will be 
occupied in an ideal GM layer at the nonequilibrium s.s. 
characterized by the global temperature f3 . Note that 
(j){e = mgz)/ N decreases monotonically with z from 1 to 
0. In fact we can write 4>(e)/ N as a cumulative probabil- 
ity 4>{e)/ N = 1 — Jq tp(e') de', and interpret the density 
■0(e) as associated to the fluidized gap. We write the s.s. 
mass profile as 



4>{z)/N 



(exp(/3mgh) — 1) 



(exp(/3mgz) + exp({3mgh) — 1) 



(12) 



To measure the profile (|12H we have implemented a sec- 
ond experiment on a GM bed (h = 21mm) with a laser 
light barrier of 10mm wide, Fig. 1(b). The voltage signal 
from the position sensitive detector runs from to 10V, 
which means a vertical window from 10mm to 0mm. We 
take measurements every 3s during 3h, where the nor- 
malized frequency count integrated from such a register 
is equal to the occupation number <f>/ N for z > fio/mg 
and to 1 — cf>/ N for z < fj,o/mg. At a fixed f e and for 
r = 10 two registers were obtained for f e = 100Hz and 
f e = 125Hz, which were integrated and normalized to get 
the corresponding profile <fi(z)/N. In Fig. 3(b) we show 
the profile and our theoretical prediction From 
this data we obtain the values (3 — 4.4 ± 0.2/iJ -1 for 
f e = 100Hz, and f3 = 10.5 ± O^J" 1 for f e = 125Hz, 
that we compare with the results of the first experiment 
Fig. 3 (a). Not only the agreement is good, but also this 
procedure allows a self-consistent test. 



Discussion. In Fig. 2 (a) we show the amplitude disper- 
sion er 2 against the velocity squared dispersion a v of the 
realizations z(t). For weak amplitude the slope a v j a z 
shows a linear behavior and the departure from a linear 
behavior is a clear evidence of the complex behavior of 
the GM bed. This indicates that for this regime it is nec- 
essary to introduce a description in terms of our theory. 

The corresponding global temperature (fcs/3) -1 for a 
fluidized gap to occur happens to be at T = 9.2 ± 0.5 x 
10 15 K, which means f e = 130Hz for T = 10, see the 
transition in Fig. 3(a). Feitosa et al. measure for a di- 
lute granular gas a range of temperatures of the order 
of T ~ 500PK (Peta Kelvin, see Fig.5 of Rcf. 0). 
This range is higher than our measurements, however 
it is in agreement with them since our corresponds to a 
weakly fluidized GM. From jjj we know that the move- 
ment of the ?7iFp when the fluidized gap appears can be 
approximated by a Brownian motion, but this descrip- 
tion changes to a more complex stochastic behavior by 
decreasing f e (for fixed T) when the temperature reaches 
T ~ 15PK. This global temperature should be under- 
stood, indeed, as equivalent to an order parameter of the 
stochastic process in the energy configuration. We re- 
mark that around the fluidization transition where the 
stochastic dynamics start to apply dz / dt occurs with 
larger rugosity. Then a proper description in terms of dif- 
ferentiable realizations z(t) is well defined in the weakly 
excited region (T > 9PK) where the dynamics start to 
be non-Markovian. For lowertemperatures (T < 9PK) 
despite of the larger rugosity of realizations z(t), the cal- 
culation of o v from a time-window of 2 second makes it 
reliable. This numerical calculation only introduces, for 
such range, a larger dispersion of the data in the curve of 
Fig. 2(b), but again in good agreement with the equation 
©■ 
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